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Abstract 

We calculated numerically the localization length of one-dimensional Anderson model with cor- 
related diagonal disorder. For zero energy point in the weak disorder limit, we showed that the 
localization length changes continuously as the correlation of the disorder increases. We found that 
higher order terms of the correlation must be included into the current perturbation result in order 
to give the correct localization length, and to connect smoothly the anomaly at zero correlation 
with the perturbation result for large correlation. 
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I. INTRODUCTION 



Electronic transport properties, the motion of electrons, in a random potential are closely 
related to the phenomenon of Anderson localization 111. The phenomena of Anderson local- 

n r n 

ization have been studied in various fields including photonics [2'], cold atoms [3], circuits |^, 
and DNA molecules (5|, |6|. Many accurate numerical approaches have been developed, by 
the quantum transfer matrix renormalization group method for finite temperature systems 
[TJ, the density matrix renormalization group method for interacting systems l8|, and the 



integral equation method for systems in the thermodynamic limit [9|, |10j, respectively. In 
this work we will study the zero energy behavior for the one-dimensional model with cor- 
related weak diagonal disorder. We first extend the numerical method we developed earlier 
in Ref. J9| for uncorrelated disorder to correlated system. Our numerical method was an 
application of the transfer matrix method ll| in localized phase in the thermodynamic limit. 
In one-dimensional Anderson model jl| with diagonal disorder is described by. 



^i_i -F Vi+i = {E - ei)i)i, 



(1) 



where hopping term is set to unity and ipi is the electron wavefunction at site-z. e^ is the 
on-site energy with a certain type of random distribution which satisfying an exponential 
correlation: {ef) = a"^ and (ejej) = cr^exp[— 12 — j|//cor] for different sites, cr^ and Icor are the 
strength and correlation length for the disordered on-site energy, respectively. Uncorrelated 
disorder is given by Icor -^ 0. Recently the anomaly around the band edge E = ±2 has 
been carefully investigated. 12| In the following we focus on the zero energy anomaly with 
exponentially correlated diagonal disorder. 

All the eigenstates are exponentially localized for one-dimensional uncorrelated disordered 
systems. 13|] The Lyapunov exponent 7 is the inverse of the localization length. It is well 
known that for the zero energy anomaly of the uncorrelated disorder system, the Lyapunov 



exponent 7 is singular at E 
clear 



a 



0. 
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15| . The physical picture behind was also 



16| . For a box distribution of uncorrelated disorder with width W and height 1/W, 



the perturbation resu 
disorder strength W 



tj;evealed that the Lyapunov exponent depends only on energy E and 



17| 



7 



96(1-^2/4)' 



(2) 



At the band center, another perturbation yielded 18|-|20| 



1 = ■ 3) 

' 105.045 ■ • • ^ ^ 

The standard variance of the disorder is a^ = W'^/12. In uncorrelated systems, order by 



order perturbation expansion in a'^ and E/a"^ has been demonstrated 10|. 

For exponentially correlated disorder, the formula for the Lyapunov exponent at finite 
energy and in the weak disorder strength limit is given by 21 1 



_ 0-2 sinhj^ 

It is straight forward to take the uncorrelated limit Icor — ;■ of formula Eq. (jlj), then obtain 

7/0"^ = 1/8 when E approaches to 0, i.e. lim lim — - = 1/8. On the other hand, if we 

stay at i^^ = 0, we should have 7/0"^ = 1/8.754 in the uncorrelated limit in accordance to 

7 
Eq. ([3]), which implies lim lim ^r = 1/8.754. Therefore, we found that the order of the 

limiting processes for E -^ and a^ — )■ can not be interchanged. It means that the point 

E = a = remains singular for perturbation expansions in o"^ and E for correlated disorders. 

The existence of strong anomalies phenomena in a correlated disorder system was pointed 



by Titov and Schomerus 22]. In this work we study the anomaly at i^^ = 0. 



II. PARAMETRIZATION METHOD 

In the transfer matrix method, Eq. ([T]) can be written as 

\ ^^ J V ^ / v^^-i 

where Tj is the transfer matrix. 

Using a parametrization method of the transfer matrix proposed in our previous work 
[9I, [lO], we will calculate the Lyapunov exponent in the thermodynamic limit within the 
localization regime. Let M^^ = TlTl_i ■ ■ ■ Ti. Then we parameterize MM* as follows 

U(0^)Mz.MiU(-^i) = ^_^^^ I , (6) 

where M* is the transpose of M and 

, cos 9r — sin Or , 

Vie,) = I . / / I • (7) 

sm Ul cos Ol 



The recursion relation of 9 in the large L limit is 

^""^""^ = i?-ej-tang, - ^'^ 

We introduce the correlations between e^ through the transformation of a group of inde- 
pendent random variables rji in terms of an identical Gaussian density distribution, 

p^{r,) = (1/v^(t) exp[-r772(T2]. (9) 

Let q = 6^^/''=°'', the exponentially correlated variable is generated implicitly by e^ = 

oo 

a/ 1 — g^ y^rji^iq'', or equivalently in the following recursive form, 

1=0 



eL = VI - q^VL + qei-i. (10) 

The three parameters X, 6, and e at a step L are what we need in order to calculate new 
parameters for the next step L + 1. 

In the localized region, the equation we obtained for the density distribution function 
p{e,e) is 

p{9,e) = ^ [ d7^de'd9'p,{r])p{9',e')S{e- ^/l^V-qe')S{-^ + ts.nd' -E + e). (11) 
sm u J tan u 

After we numerically solve this equation, the Lyapunov exponent 7 can be calculated through 
the following formula, 

7 = - / dr]dedep^{7])p{e, e) ln[l - (E - ^/l - q^r] - qe) sin 20 + {E - ^Jl^i] - qef cos^ e\. 

(12) 
By defining the distribution function p{Q\ which is similar to the one in the uncorrelated 
diagonal disorder case, 

p{Q)= fp{9,e)de, (13) 



we obtain the same simple relationship between p{6) and the Lyapunov exponent, 

7 = - I p{e)\n\tan9\d9. (14) 



If we take the limit of Icor — ;■ in the present correlated disorder situation, the equations 
of p{6) in the uncorrelated disorder case will be recovered 9|, [lO|. However, p{9,e) is not 
exactly the product of p(6')p^(e) in this limit. 




FIG. 1: Distribution p{9, e) for E = Icor = o" = 1. The forty lines in direction are evenly spaced 
in the region [— 7r/2,7r/2]. The forty lines in e direction are scaled to display a better global view. 

We use the Gaussian distribution p^ to solve Eq. fITT]) and to calculate 7 numerically. 
This method is very efficient to yield high precision results for various disorder correlation 
length Icori disorder strength o"^, and energy E in the thermodynamic limit. Fig. [T] and 
Fig. [2] are shown the calculated distribution functions p{6, e) and p{6), respectively. In these 
calculations we have set a relative precision 10"^'' for p{6,e). Similar distributions were 
calculated recently [23| in the dichotomous correlated disorder case. 

The structure of the joint density distribution of 6 and e is demonstrated by the p{6, e) 
of ii^ = Icor = cr = 1 in Fig. [TJ The distribution is not so complicated to perceive, but it 
can not be decomposed into a direct product of a density distribution for 6 and a density 
distribution for e. Two curves for p{6) oi E = 1 are given in Fig. |2l One of the curve 
with Icor = o" = 0.01 has very small disorder strength a and very small disorder correlation 



Icor- This curve can be approximated very well by the expression 



'or distribution p{9) of 



,u 



uncorrelated disorder at a finite E in the weak disorder limit 10|, llSl, l20|] , 

,(e) = —^^^l^—„ (15) 

7r(l — cos/ism/t7j 
where cos/x = E/2. Another curve with Icor = o" = 1 is not in the case for small disorder 
strength or small disorder correlation, which is different from the curve in small disorder 
strength or small disorder correlation. 

The Lyapunov exponent 7 (/cor, cr) is then calculated by using the two curves shown in 
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FIG. 2: Distributions p{9) at E = 1. The full line is for l^or = cr 

hor = o" = 0.01. Each point on the dotted line is differed from the analytical curve 

within no more than a relative error 10~^. 



1; and dotted line for 
7r(2-sin26') 



Fig. Elat E = 1. We obtain 7(1, 1) = 0.1252 and 7(0.01, 0.01) = 0.00001667. The direct cal- 
culated results from formula Eq. (g]) yield 7(1, 1) = 0.09587 and 7(0.01, 0.01) = 0.00001667. 
We see the numerically calculated Lyapunov exponent for the finite energy and in the weak 
disorder strength limit is well predicted by formula Eq. (jlj). 

The situation for zero energy is different compared to that for the finite energy. There is 
no analytical result obtained so far for the zero energy anomaly in the presence of correlated 
disorder in the weak disorder limit; nor the formula predicting the Lyapunov exponent for 
a finite correlation length. Our method is a good choice to perform calculation in these 
situations. 

III. ANOMALY AT £; = 



We will investigate how the localization length changes as the correlation of the disorder 
varies at i^^ = in the weak disorder limit. At finite energy, the disorder strength a and 
the correlation Icor are decoupled in function 7(/cor,o") in Eq. (^. Since Eq. (jlj) was derived 
without any limitation on the magnitude of the correlation length, with the help of the 




0.34 



FIG. 3: Distribution p{6,e) for E = 0, tanh 2] — = 0.475, and a = 0.1. The forty lines in 9 
direction are evenly spaced in between [— 7r/2, vr/2]. The forty lines in e direction are scaled to give 
a better global view. 

Lyapunov exponent 7(0, a) for uncorrelated disorder, the ratio ■y{lcor, c")/7(0, cr) = tanh ^f— 
might be exactly held for any Icor- At E = anomaly, even if we keep only a^ term 
in the weak disorder limit, it is not known whether higher order terms from correlation 
exists beyond the perturbation result. To answer this question, we compare the numerically 
calculated result with the perturbation one given by Eq. i^ at E = 0: 

^2 1 



a 



Iv 



— tanh 



2L 



(16) 



The small quantity related to correlation in 7p can be considered in two limit cases. In the 
short correlation length limit Icor -^ 0, i.e. tanh jt^^ > 1, the small quantity for expansion 



is 1 — tanh 



^I'cor 



2e ^/'<=°'-; whereas in the large correlation length limit Icor -^ 00, the small 
quantity for expansion is tanh 757^ itself, tanh ^^y^ ~ iw~- Therefore, we will calculate for 
a group of different correlations with tanh 1^^ close to zero as well as to one. In order to 
neglect the contribution from the higher order terms of a in our calculation in the weak 
disorder limit, we will calculate only for small a. It is sufficient to keep three significant 
digits for the Lyapunov exponent '^{Icor,'^)- 

To demonstrate the anomalous behavior at -E = 0, we plot in Fig. [3] the distribution 
p{0,e) for i? = 0, tanh 2f— = 0.475, and a = 0.1; and in Fig. IHthe distributions p{6) for 




.925 



FIG. 4: Distributions p{9) for E = and a = 0.1. The forty lines in 9 direction are evenly spaced 

0.025, 0.075, 



in between [— 7r/2,7r/2]. The twenty lines in x direction are for x = 1 — tanh 2j^ 
0.125, . . ., 0.925, 0.975, respectively. 



E = and o" = 0.1 with x 



tanh 



^/•cor 



0.025, 0.075, 0.125, . . ., 0.925, and 0.975, 



respectively. It shows clearly in Fig. [3] that the joint distribution for 6 and e has some inner 
structure. We have observed the flattening of the distribution p{6, e) when increasing the 
correlation length l^or in the weak disorder limit. The flattening will not be presented in 
Fig. El 

The flattening of p{6) can be seen in Fig. HI In the flgure, when correlation length is 
small, the distribution p{6) turns out to be similar to the distribution for the uncorrelated 
disorder|lo|,y,y: 

where K is the complete elliptic integral of the flrst kind. As the correlation increases, 
we see that the distribution p{6) flattened towards l/vr. Let's take x = 1 — tanh 



^Icor 



as 



a new parameter of the correlation in disorder. In the limit of x — ?■ 0, which corresponds 



to the uncorrelated limit l^or -^ 0, both the anomalous distribution l/\/3 + cos 4^ and the 
anomalous Lyapunov exponent 7 = cr^/8.754 of the uncorrelated disorder will be recovered. 
In the limit of x — t- 1, which is equivalent to the large correlation limit Icor "^ 00, p{9) = l/vr 
will correctly give a zero Lyapunov exponent. 
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FIG. 5: The Lyapunov exponent 7 for E = and a = 0.1. The variable x used for different 
correlations is x = 1 — tanh tt^— • The function is w = 7/(0"^ tanh ttt^)- When x is close to zero, y 

Zlcor ZLcor 

is close to 1/8.754; and when x is close to one, y is close to 1/8. 

IV. HIGH ORDER TERMS OF CORRELATION AT E = 



Now we analyze the contribution from higher order terms of the correlation in the weak 
disorder limit. In Ref. 25| the authors gave analyses, which cover not only the localization 
length, but also all the higher moments of the distribution of the Lyapunov exponent for 
uncorrelated finite systems. For correlated systems we expect the deviation from Eq. f fT6|) 
comes from higher order terms of correlation too. 

In the perturbation result 'y„ in Eq. f fT6l) . by using variable x = 1 — tanh ^rj^ to denote 
the correlation, we see that 7^ included the first order correction of small x when x — )■ 
and also the first order correction of small 1 — x when x — ;■ 1. 7^ has included only the 
first order term. From the discussion on i? = anomaly in the previous section we know 
that 7/(cr^tanh 2|^) = 1/8.754 for x — )■ 0, while 7/((T^tanh gf— ) = 1/8 is predicted by 
perturbation result for E ^ 0. The question on how 7 really behaves at ii^ = is still 
not answered: whether 7/(0"^ tanh tjt^) = 1/8.754 always holds, or there is a crossing to 
7/(0"^ tanh 2f—) = 1/8 as Icor increases. We plot Fig. E] to answer this question. 

In Fig. El the Lyapunov exponent 7 for E = Q and cr = 0.1 are presented. We plot for dif- 



ferent correlations by using the parameter x = 1— tanh ^i^, and we plot y = 7/(0"^ tanh 2r~) 
as the function of x. When x is close to zero, y is close to 1/8.754; and when x is close 
to one, y is close to 1/8. We see a crossover between the anomalous value 1/8.754 and 
the perturbation result 1/8. In the weak disorder limit, besides the term tanh2f— , there 
are higher order terms in x or 1 — a; from the correlation. The higher order terms connect 
smoothly the anomalous 1/8.754 at zero correlation with the perturbation result 1/8 for 
large correlation length. 

The physical picture is rich behind a finite magnitude of a and a large correlation length. 
The a in Fig. \5\ is not a small enough disorder strength. The higher order terms in o"^ 
contributes when x approaches one in Fig. O We have calculated for much smaller a and 
confirmed that the contribution of higher order terms in a^ goes to zero in the weak disorder 
limit. Our observation suggests further perturbation investigations. 

To numerically provide the next leading term of the correlation closed to the uncorrelated 
limit, we fit 7 for x close to zero in Fig. [6l In Fig. |6] the Lyapunov exponent 7 for E = and 
a = 0.01 is plotted, y represents the difference between the Lyapunov exponent for a finite 
correlation and for zero correlation: y = 7 (/cor, o")/(o"^ tanh ^j^) — 1/8.754. The variable x 
used for different correlations is a; = 1 — tanh 777^. We obtain a fitting line y = 0.01533a;. 
Therefore the perturbation expansion of 7 to the sub-leading order of the correlation in 
power of X is obtained, 

.= (1 + 0.1342.)^ tanh J-. (18) 

In Fig. [5] it is clear that higher order terms contributes when x is even bigger. In the weak 
disorder limit, when 1 — x close zero, the next order term in the correction factor used to 
multiply to 7^ in Eq. flT6|) is (1 — x)^. 

V. CONCLUSION 

In summary, we calculated the inverse localization length in one- dimensional Anderson 
model with correlated diagonal disorder. We obtained numerically the curve of the inverse 
localization length for correlations at zero energy in the case of weak disorder. A nonsingular 
curve was obtained for different correlation lengths in the weak disorder limit at zero energy. 

The variable used to plot the unifying curve is tanh t^t—, which has correspondence to the 
Poisson process of the phase accumulation. The inverse localization length will be singular 
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FIG. 6: Fitting of the Lyapunov exponent 7 for E = and a = 0.01. y is certain the 
difference between the Lyapunov exponent for finite correlation and zero correlation: y = 

iQ-cor: (T)/c'"^/tanh 2] 1/8.754. The variable x used for different correlations is x = 1 — tanh 2j — • 

The fitting line is y = 0.01533x. 

as the function of other variables as Icori ^/hor, or e"^/'""*". We suggest further studies on 
the inverse locahzation length in perturbation expansions or functional expansions with the 
parameter tanh 2f— • We have obtained numerically in this work the next leading term for 
comparison. 

We also saw rich behavior for finite disorder strength and large correlation length. A 
unifying description of the band center anomaly and the correlated disorder will be very 
interesting. 
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